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We extend our Pade-aided analysis of the nonperturbative renormalization of nucleon-nucleon 
scattering to the case of coupled channels. 
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Since the significant effective field theory (EFT) approach to nucleon systems [HQ, there has been creating 
controversies about the consistent renormalization and EFT power counting in nonperturbative regime 0,0). 
6JT)' The main difficulty is due to the nontrivial prescription dependence developed in nonperturbative regime, 

where some wisdoms about renormalization established within perturbative regimes cease to apply directly 
Therefore, the nonperturbative prescription dependence must be removed through imposing appropriate 
boundary conditions which are usually implemented through various forms of data fitting either explicitly or 
implicitlvfq. |7L fsL [oL Hol. HH Ha . Il3| ■ In more recent literature, nonperturbative counter terms, or equivalently, 
nonperturbative parametrization of the renormalization prescription has been the main focus|14l. |15|. That 
is, due to the difficulty in treating the issues in nonperturbative regimes, the key issue is to find more efficient 
parametrization of nonperturbative prescription. 

In this regard, we have performed an analysis and treatment of the nonperturbative prescription depen- 
dence basing on Pade approximant of an important nonperturbative factor of T-matrix for nucleon-nucleon 
(NN) scattering in uncoupled channles [Til [ttI fill]. In this report, we will extend our analysis to coupled 
channels. Some general theoretical and technical issues associated with coupled channels will be addressed 
first, then we will illustrate our method in 3 r>3— 3 Gs. 

The object under consideration is the T-matrix for nucleon-nucleon scattering processes at low energies. 
According to Weinberg's proposal[l^|, this T-matrix should be solved from Lippmann-Schwinger (LS) equa- 
tion with the potential to be systematically constructed using xPT as the low energy effective theory of 
QCD. For coupled channels, such LSE's read, 
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Go{k;E+) = E+ _ k2/M , E+^E + l e, (1) 

with E being the nucleon energy in the center mass frame, M the nucleon mass, p' — |p'|, p — |p|. The 
bold-faced capital letters represent the 2x2 matrix-valued objects in the angular quantum number space. 
The convolution is understood as already regularized and/or renormalized in an unspecified prescription in 
order to make our discussions generally valid. Following Refs.[l^,[l3,[lB|, the above LSE for coupled channels 
could be transformed into a compact and hence nonperturbative parametrization of T-matrix as below, 

T- V,?>; E) = V-^p^p; E) - g(p',p; E), (2) 



(p',p; E) = V(p',p; E) + J G (fc; E+)V(p', k; E) x T(k,p; E), 



g{p',p;E) = V- 1 {p',p ] E) x 



J G (fc; E+)V(p', k- E) x T(k,p; E) 



y.T- l {p',p-E) (3) 



where the factor Q assumes all the 'loop' processes generated by V in the field-theoretical terminology. 
Making use of the X-matrix formalism, the unitarity of such compact T-matrices follows immediately (lg, 
UtI . [l8j . It is also easy to verify the inverse relation in the coupled channels: 

Txr'=|vi J G V x T^j x T" 1 = V x T" 1 + V x Q = V x (T^ 1 +g)=V x V" 1 = I, (4) 

T" 1 x T = (V" 1 - Q) X T) = V" 1 x T - g x T = V" 1 x — J G V x T^j = V" 1 x V = I, (5) 



* corresponding author. 
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with I denoting the 2x2 unit matrix. Obviously, in terms of V and Q, T are nonperturbative objects. 
It is interesting to note that in terms of the standard parametrization of S'-matrix, 



cos 2t (p) exp[2i*y j (p)] i sin 2t 3 (p) exp [i{5) 3 _ 1 (p) + 5^ (p))] 

ism2e,(p) exp[i(S] 3 +1 (p) + ^(p))] cos2e,(p) exp [20]^^)] 



i — -T, 

2tt ' 



(6) 



we could find that, 
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(7) 



Thus, the imaginary part of T is simple and proportional to a unit matrix, the off-diagonal entries of T 
are real numbers only. The unitarity now reads: T _1 — (T -1 )* = i^-I- Only the real part of T" 1 (or Q) is 
subject to nonperturbative renormalization and hence the unitarity is not affected by renormalization. The 
prescription dependence is exclusively contained in the real part of the factor Q . 

The motivations and plausibility of employing Pade approximant to Q were already demonstrated in 
Refs.[ll,[l2l 

[l8j . In this approximation, the unitarity of T-matrices is automatically preserved, a virtue that 
is welcome in hadron physics (20|. Another important virtue is its generality and flexibility in parametrizing 
the prescription dependence, avoiding being stuck in or confined to a special prescription that might not be 
quite compatible with physical boundaries. 

Now it is clear that Q is a 2 x 2 matrix with the diagonal entries being complex, while the off-diagonal 
ones real. The Pade approximant is applied to the real part of any matrix element, diagonal or off-diagonal 
(p = VME): 



Re{£(p)i+A,i+A'}|pade 



J2k N l+A,l+A';kP 2k 
Efe' A+A,i+A';fc'P 2fc ' 



A, A' = 0,2. 



Evidently, in Pade approximant, the prescription dependence is contained the parameters [N...-k, D...-k'] 
other words, [N...-k, D. ..■/.'] serve as approximate parametrization of the nonperturbative prescription. 
Then the renormalized T could be approximately parametrized in nonperturbative regime as follows, 



(8) 
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a (p; [g..., C...}; [N..,D..]) = V-\p; [<?.., C...]) - 



Efc Ni+AJ+A'-MP 
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Z)fe' D l+A,l+A';k'P 
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Simple and coarse as it is, such nonperturbative parametrization of the T-matrix contains all contributing pa- 
rameters: EFT couplings [<?..., C...], and renormalization prescription parameterized in terms of [N...-k, D. ..;/->]■ 
As a byproduct, the Pade parameters allows us in principle to effectively imitate any renormalization pre- 
scription of T through corresponding definition of [N...-k, D...-k']- 

For our approximation to be sensible, the Pade parameters must be appropriately determined. Then their 
magnitude orders should be in accordance with EFT power counting. For a general EFT power counting, 
one may expect that: 



N.. 



D. 



N... i0 ' D....o 



(io) 



with A(~ 500MeV) being the upper EFT scale, fj, the typical EFT scale, here, say, ~ (10, 100)MeV. F, f 
are some counting functions. Note that the reflection of EFT power counting in the factor Q is completely 
nonperturbative, in sheer contrast to the conventional understandings. 

Occasional large deviation from such rules should be due to unnatural behaviors of the NN scattering. 
The EFT approach would be indeed problematic only if no power counting scheme could be sensibly realized 
in any renormalization prescription. In other words, the failure of some power counting schemes does not 
imply the very failure of the EFT approach. 

In general, the intrinsic scales involved in the Pade parameters should be A and p.: 



N...., 
N..., 



D. 

w. 



//"A 1 



a e (0,1.0), fc> 0. 



(11) 
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As dim[Cf] = 2, we choose dim[D... ; o] = 2 and hence |Z?... ; o|^ ~ /i Q A 1_Q , a £ (0,1.0). This is because Pade 
parameters are in fact functions of both EFT couplings and renormalization scales or constants. 

Now, to explore physics using the parametrization given in Eq.((5]), we must fix the prescription or Pade 
parameters through imposing appropriate boundary conditions. To this end, as in most literature, we fit to 
the PWA[2~fl| data for the phase shifts and mixing angle in the low energy ends, say the kinetic laboratory 
energy TLab(= 2E) £ (0, 50)MeV. In coupled channels, one must fit three sets of Pade parameters for the 
phase shifts and mixing angle at the same time. In order to see the main points or rationalities of the 
Pade-aided analysis of coupled channels, wc will work with the simplest and hence coarsest cases of Pade 
approximant, i.e., the constant Q factor, 



Alternatively, with such choice of pade approximant, we wish to probe the most important scales in the 
nonperturbative factors G, in order to see if there would be significant deviation from the EFT power 
counting, or abnormal numbers. 

In this short report, we pick up the 3 D^— channels for illustration where up to next-to-next-to-lcading 
order the potentials contain no extra contact terms to be determined first. This is also the highest coupled 
channels where at least one channel, 3 £>3, is not perturbative, which could be seen below. As before, we 
employ the potentials and couplings given by EGMJt]]. One could well employ other sets of definitions for 
comparison. Such works will be carried out in the future. We also note that we deliberately work with low 
precision in order to save computer workloads: round up to the first two digits. 

The numerical results are summarized and presented in TABLE I and Fig[T] The nonperturbative renor- 
malization does significantly improve the phase shift predictions for the 3 Z?3 channel in comparison with 
the inferior perturbative ones as depicted in Fig[2j where the perturbative predictions even produce wrong 
sign of the phase shifts. For Ssq 3 we find similar but less significant improvement for lab energies below 
lOOMeV, as the perturbative predictions for all the three orders already deviate from the PWA data from, 
say TLabsjm,75MeV, though less significant. Such simple results already means that only the prescription 
fixed through physical boundaries could reliably describe physics. (One could try other rather different 
values of g .. and see that the phase shifts and mixing angle thus obtained are nonsense.) As the p or en- 
ergy dependence in Re(G) is totally discarded here, the predictions are doomed to fail as E is higher, say, 
TLab > lOOMeV, which is evident from Fig[T] and FigO Also the trend that the predictions improve order 
by order is not clear here. To see this trend, more sophisticated Pade approximation and hence heavier 
workloads are required, which will improve the predictions in many respects. Further works along such lines 
are in progress and will be reported in the near future. Here, we are merely content with illustrating the 
plausibility of Pade approximant to the factor Re(C/). Although the Pade approximant adopted here is very 
coarse, the main virtues in using such relatively more analytical and controllable approach are still quite 
significant from the simple numerical analysis given in the figures and tables. 

We still need to show that the Pade parameters obtained via fitting, here, g .., follow the rules described 
above in Eq. fTTj) or (fTU)l . To this end, we have computed square roots of the absolute values of g .. and 
listed them in Table I. From Table I one could see that the scale extracted from the coarse nonperturbative 
approximation lies between 10 and 200MeV, just in the range described by Eq. tTiTj) , that is, (10, 500)MeV. 
In terms of a, we have a £ (0.26,0.91) C (0, 1.0). In other words, through the Pade approximant of the 
factor Q, the scales involved in the nonperturbatively renormalized T- matrices do not fall outside of the 
EFT's scope. 

At this stage, we may conclude that the Pade-aided approximation to renormalized T-matrices also works 
in coupled channels. 

In summary, we extended our Pade-aided analysis into coupled channels. Primary numerical analysis 
showed that such treatment also works in the coupled channels. The results also exhibit that intensive and 
extensive studies are needed to further develop this promising treatment for investigating various issues in 
the EFT approach to nucleon-nucleon systems. 
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TABLE I: The simplest Pade parameters for Re(C?) (/MeV 2 ) fitted at different chiral orders. 
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FIG. 1: Predictions of S3 D3 , S3 Ga and 63 versus lab energy T^b in MeV with the simplest Pade, with solid line for 
PWA, dotted lines for LO, dashed line for NLO and dot-dashed line for NNLO. 
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FIG. 2: Perturbative predictions for S3 Da , S 3 g 3 an d £3 versus lab energy TLab in MeV, i.e., G=0. Conventions are as 
in FigfT] 



